Theory of Kinetic Partitioning in Protein Folding 
(With Application to Prions) 

V. I. Abkevich, A. M. Gutin, and E. I. Shakhnovich 

Harvard University, Department of Chemistry 
12 Oxford Street, Cambridge MA 02138 
(February 1, 2008) 



Correspondence to E. Shakhnovich 
Dept of Chemistry Harvard University 
12 Oxford Street Cambridge, MA 02138 
tel 617-495-4130 
FAX 617-496-5948 
Int ernet : eugene @ diamond . harvard . edu 

Keywords: prion protein; conformational transition; protein folding 



1 



Abstract 

In this paper we study the phenomenon of kinetic partitioning when a 
polypeptide chain has two ground state conformations one of which is more ki- 
netically reachable than the other. This question is relevant to understand the 
phenomenology of prions, proteins which exist in the cell in non-pathogenic 
a-helical conformation but under certain circumstances may transform into 
pathogenic PrP Sc state featuring increased /3-sheet content. We designed se- 
quences for lattice model proteins having two different conformations of equal 
energy corresponding to the global energy minimum. Folding simulations 
revealed that one of these conformations was indeed much more kinetically 
accessible than the other. We found that the number and strength of local 
contacts in the ground state conformation is the major factor which deter- 
mines which conformation is reached faster: The greater the number of local 
contacts the more kinetically reachable a conformation is. We present sim- 
ple statistical-mechanical arguments to explain these findings. The presented 
results are in a clear agreement with experimental data on prions and other 
proteins exhibiting kinetic partitioning. 
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INTRODUCTION 



The study of prion biology and diseases is a new emerging area of biomedical investi- 
gation (reviewed by Prusiner, 1992). Prions are proteinaceous infectious particles that are 
composed largely, if not entirely, of an abnormal form of the prion protein (PrP) desig- 
nated, in the case of scrapie, PrP 5c (Prusiner, 1991). The term "prion" was introduced by 
Prusiner (1982, 1991), who has shown that prions are unique among all infectious pathogen- 
esis may be both inherited and transmissible. Six diseases of animals and four of humans 
are caused by prions (Prusiner, 1992). The initial hypothesis was that PrP Sc was derived 
from normal cellular prion protein (PrP c ) by a post-translational process(Borchelt et al, 
1990, 1992; Caughey & Raymond, 1991; Taraboulos et al, 1992). Attempts to identify a 
post-translational chemical modification have been unsuccessful (Stahl et al, 1993). Recent 
structural studies demonstrate that PrP 5c and PrP c have the same chemical structure but 
dramatically different conformations (Pan et al, 1993). This suggests that prions have two 
low free energy states in both of which they can be stable during the lifetime of a protein, 
with pathogenic form PrP 5c being much less soluble than than cellular PrP c . The prion 
puzzle has two aspects. First, it is the reason and the character of transition from cellular 
to the pathogenic form. The most intriguing aspect of it is unusual infectivity by PrP Sc 
which clearly suggests that intermolecular interactions are likely to play an important role 
here. A possible mechanism of transition from PrP c to PrP 5c has been suggested recently 
(Cohen et al, 1994) and its implications for prion pathogenesis were discussed. 

The focus of this paper, however, is on another important aspect of prion puzzle, namely: 
Why after (or concomitantly with) the synthesis of their primary structure, prions do not 
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fold into the conformation of PrP , at the first hand? 

One (trivial) explanation may be that PrP*" is much more thermodynamically stable 
than monomeric PrP 5c . However, the experimental data do not seem to support this view. 
One (indirect) evidence is that PrP c is much more susceptible to proteolysis than PrP 5c . 
Pan et al (1993) suggested that PrP * is a kinetic trap. 

If PrP c and monomeric PrP 5c have comparable thermodynamic stabilities and the par- 
titioning into normal and pathogenic conformations was driven by thermodynamic rules, the 
yield of pathogenic species after synthesis of primary structure should be large. However, 
this is not so, and dominant form in normal, uninfected cells is PrP c . Apparently, this 
may be so because PrP c conformation is more reachable kinetically. This links prion phe- 
nomenology to the most fundamental aspects of protein folding since they provide a clear 
(but not unique, see below) example of kinetic partitioning when a conformation becomes 
dominantly populated for kinetic, rather than thermodynamic, reasons. 

In this paper, we, intrigued by the mystery of prions, study the phenomenon of kinetic 
partitioning in protein folding. Currently, folding simulations are feasible only in the realm 
of simplified lattice and off-lattice models, which already provided useful, experimentally 
testable, insights into such general principles of folding as nucleation mechanism (Abke- 
vich et al, 1994b; Shakhnovich et al, 1996; Fersht, 1995; Itzhaki et al, 1995) and folding 
intermediates (Gutin et al, 1995a; Mirny et al, 1996; Fersht, 1995). 

Our approach to study kinetic partitioning, relevant to prion folding, aims to simulate 
it in a simple, yet, nontrivial, model and seek full understanding of this phenomenon within 
the framework of the model. Then using the model results as an example where kinetic 
partitioning is understood, one can apply this as a lead to thinking about real systems and 
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as a tool to understand existing and plan future experiments. Of course such approach may 
be useful only if there is a generic reason for kinetic partitioning, which may be applicable 
not only to prions, but to a wider range of systems. In this paper we reveal a possible 
physical mechanism of kinetic partitioning which is due to global structural properties of 
native conformations and we argue that it may explain certain observed features of prion 
folding as well as other proteins where kinetic partitioning is likely to be important. 

In our model study we would like to focus on the kinetic aspect of prion folding. There- 
fore, while it is not known which conformation, PrP c or PrP 5c , has lower free energy, in the 
model, we assume that two "native" conformations have equal (free) energies and seek the 
reason why one of them is more reachable than the other. 

Specifically, in this work we design sequences having two different conformations of global 
energy minimum and study their folding. Our results provide structural clues to kinetic 
partitioning explaining, why in a system with two ground state conformations, one of them 
can be more kinetically reachable than the other. 

We compare the results of our study with phenomenology of prions and other proteins 
having kinetic partitioning. 

THE MODEL 

We represent proteins as self-avoiding chains on an infinite cubic lattice, such that co- 
valently linked residues occupy neighboring lattice sites. The energy of a conformation is 
the sum of energies of pairwise contacts between monomers. Two monomers are defined to 
be in a contact if they are neighbors on the lattice and not connected by a covalent bond. 
The energy of a contact depends only on the identity of the two amino acids involved. The 
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interaction energies for amino acid pairs are determined from the statistical distribution of 
contacts in real proteins (Miyazawa & Jernigan, 1985, Table VI). 

Our approach to folding simulations requires first to choose a target conformation to be 
the native and design amino acid sequences that fold to, and are stable in, this conformation. 
It was shown previously (Goldstein et al, 1992; Shakhnovich & Gutin, 1993; Gutin et al, 
1995b) that such a sequence should render the target native conformation as a pronounced 
global energy minimum. This is the criterion used in our design algorithm (a Monte Carlo 
optimization in sequence space). The details of this algorithm have been published elsewhere 
(Abkevich et al, 1995). The quantity that is minimized in this design is relative value of 
native energy comparatively to the non-native conformations (Bowie et al, 1991; Gutin et 
al, 1995b): 

^ E na t E av 

a 

where E nat is energy of the native conformation, E av is the average energy of compact 
non-native conformations with corresponding dispersion a. 

Our aim in this work is to design sequences that exhibit prion-like behavior, i.e. which 
have more than one conformation as a global energy minimum. Correspondingly, the design 
procedure should be modified to achieve that goal. To this end we choose two different 
conformations and design such amino acid sequences that have low energy in both confor- 
mations. The essential parameter which is minimized in this design is 

Z = Z 1 + Z 2 + (Z 1 - Z 2 f (2) 

where Z\ and Z 2 are relative energies of the model polypeptide in the first and second target 
conformations, respectfully. The last term in the equation (2) is introduced to ensure that 
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designed sequence have close energies in target conformations. 

RESULTS 

The first model protein we studied is a 36-mer chain which has equally low energy in two 
target conformations (Fig.l) arbitrarily chosen out of more than 84 million possible fully 
compact conformations (Pande et ai, 1994). The sequence was designed to have low energy Fig.l 
in both of these conformations (Fig.lc)[|. 

In order to simulate folding, we used the standard Monte Carlo method (Hilhorst & 
Deutch, 1975). Different simulation runs begin from different random coil conformations. 
In the process of simulation designed protein always reached each of the conformations shown 
in Fig.l. Further, in a long Monte Carlo folding run (10 9 steps) we did not observe any con- 
formations having energy lower than these two which is indicative that these conformations 
are the lowest in energy. 

The mean first passage time (MFPT) | into the structure shown in Fig. lb is found at 
least 20 times longer than into the one shown Fig. la in a wide range of temperatures. To 
verify this result we measured MFPT for ten nonhomologous designed sequences that have 
close low energies in both conformations shown in Fig.l and we consistently observed much 

x One may notice that sequences shown on Fig.lc and Fig. 2c have unusually large content of W 
and M. This is due to the fact that our design program does not take into account geometrical 
properties of amino acids and their natural occurrence but only energy of interactions between 
them, which is determined only approximately. The strongest hydrophobic interactions in the 
energy set used in this work (Miyazawa & Jernigan, 1985, Table VI) is a one between tryptophan 
and methionine. That is why these amino acids dominate in our designed sequences. 

2 MFPT in all experiments reported in this paper was calculated by averaging over 50 folding runs 
each starting from different random conformations. 
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faster folding into the conformation shown in Fig. la. This result was unexpected, and the 
first guess was that one of the conformations shown in Fig.l is more kinetically reachable 
for geometric, or algorithmic reasons. To test this hypothesis we designed sequences were 
designed having each of the structures shown on Fig.l as their unique native conformation. 
MFPT for these sequences at their folding transition temperature into their respective native 
conformations were approximately equal (data not shown). This ruled out the simple, but 
mostly artificial, possibility mentioned above. 

Difference in folding time into conformations with completely different structure but close 
energy was also observed for random model polypeptides. Studying folding of ten random 
27-mers with the same amino acid composition we found that nine of them fold into their 
lowest energy conformation in about 10 6 Monte Carlo steps. However, the tenth sequence 
folded into the conformation with the lowest energy in more then 10 s steps. However folding 
time at the same temperature into a completely different conformation with only slightly 
higher energy is 20 times faster. Again such a difference in folding times appeared to be not 
just due to geometrical inaccessibility of the conformation with the lowest energy. Folding 
of the sequence which was designed to have this conformation as the global energy minimum 
was fast. 

It seems that when random or specially designed sequence have close energy in two 
conformations folding time into these conformations is often quite different. This can be 
similar to what is found for prions which fold into their native conformation and stay in 
it during life time of a protein being digested by a protease before reaching the alternative 
stable conformation. 

How can the existence of one low energy state influence folding into another? The 
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assertion that one conformation can play a role of a kinetical trap (Abkevich et al, 1994a) 
is not sufficient to explain the observed behavior. At studied temperatures random 27-mer 
is not stable in its low energy conformation, and polypeptide spends most of its time in the 
unfolded state. Since the difference in folding times remains even when the ground state is 
unstable one should seek the reason of this difference in the propoerties of the unfolded state. 
It is known that in unfolded state strong local contacts prevail (DeGennes, 1979; Grosberg & 
Khohlov, 1994). The ground state conformations are the lowest in energy and thus ground 
state contacts are the strongest on average. So among contacts even in a unfolded state those 
which are the same as in the ground state, especially local ones, will dominate. Comparing 
the structure shown in Fig. lb with the one shown in Fig. la we found significant difference 
in the number of local contacts (contacts between monomers i and i + 3). There are six local 
contacts in the structure shown on Fig. lb and eleven such contacts in the structure shown 
on Fig. la. 

When we compared the lowest energy state of a random 27-mer, for which kinetic par- 
titioning was observed, with the alternative low-energy conformation, which was reached 
much faster, we found that both conformations have the same number of local contacts. 
However, average energy of local contacts was lower in the conformation faster reachable 
conformation (-0.27) than in the native one (-0.17). 

Thus one can suggest that the number and strength of local contacts are the major 
factors determining which conformation will be reached first. To pursue this lead further, 
we chose another couple of structures for design: the one with many (seventeen) local 
contacts (Fig. 2a) and another one without local contacts (Fig. 2b). A sequence designed to Fig 
have equally low energy in these structures is shown on Fig.2c.If our hypothesis is correct, 



this sequence will fold into the structure shown on Fig. 2a faster than into the one shown 
on Fig.2b, and the ratio of MFPT should be greater than for structures shown on Fig.l. 
The results of simulations fully confirm this prediction: At temperature when folding is the 
fastest, the MFPT from random coil into the conformation shown on Fig. 2b is equal to 
5.6 • 10 6 Monte Carlo steps, whereas MFPT into the conformation shown on Fig. 2b is equal 
to 7.8 • 10 8 : more than 100 times longer. 

To show that this result is generic and does not depend on the parameter set, we extended 
analysis using so called Go model (Taketomi et al, 1975) in which all native interactions 
are set equal and attractive and all other interactions are set to zero. Apparently, such 
model does not consider any interactions except the native ones, and hence it is somewhat 
unphysical. However the Go model has some important advantages for the problem under 
investigation. Firstly, all contacts in both ground states have the same energy, and thus total 
energy of ground state contacts is simply proportional to their number. Secondly, many 
factors which influence folding rate in more realistic sequence models (dispersion of energy 
of native contacts (Abkevich et al, 1996), stability of the specific nucleus (Abkevich et al, 
1994b; Fersht, 1995) and so on) are not important for the Go model. Further, the MFPT at 
the optimum temperature for the Go model is substantially less than for a sequence model. 
This allows one to overcome considerable computational difficulties of MFPT determination. 

We compared MFPT for the Go model into the conformations shown in Fig. 2 at different 
temperatures (Fig.3). We found that folding into the structure shown in Fig. 2a is at least Fig. 3 
100 times faster than into the one shown on Fig. 2b. 

The analysis of a typical folding trajectory for a long Monte Carlo run at folding transition 
temperature (Fig.4) provides further insight into the origin of such a pronounced difference 

10 



in folding rates into two different ground state structures. Fig. 4 

For convenience of subsequent discussion we denote the structure with many local con- 
tacts (Fig. la, Fig.2a) as Ni and the structure with a few local contacts (Fig. lb, Fig.2b) as 
N 2 . It can be seen from Fig.4 that properties of unfolded state give rise to the differences in 
folding into Aq and N 2 . The chain rapidly transforms into the set of (unfolded) states with 
significant structural similarity with Aq, having about 50% of the contacts in common with 
this ground state conformation. In contrast, the unfolded conformations bear very little 
structural similarity with N 2 and only after almost 3- 10 7 MC steps it reaches conformations 
with sufficient number (and location) of N 2 contacts to enable nucleation and subsequent 
rapid folding into that ground state conformation. This is due to the fact that local contacts, 
numerous in Aq, are more favorable than other possible contacts. After these contacts are 
formed, structural similarity with Aq is enforced. This facilitates faster folding into Aq. In 
order to nucleate folding into N 2 , the most stable local contacts in the unfolded state must 
be broken. This process is uphill in free energy and therefore requires longer time. 
This qualitative picture can be represented by the following scheme: 

t 

U 2 ^N 2 

Here U\ and U 2 - unfolded states in which dominate contacts from the corresponding 
ground states. U\ might be similar to PrP* state (Cohen et al, 1994). As was noted earlier, 
folding from the state Lq into the state N 2 firstly requires to tear contacts in U 1 which are 
common with Aq but not with N 2 . This corresponds to the transition from U\ into U 2 . 
The same is also true for the states U 2 and Aq. The essential features of the free energy 
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landscape of prions are summarized on a schematic diagram shown on Fig. 5. Fig. 5 

Even if the rate of folding from U\ to N x is equal to such from U 2 to N 2 ^, the folding time 
from the random coil state into states N\ and N 2 will be different if equilibrium constant 
K = [Ui]/[Uz\ between states U\ and U 2 is not unity. If K <C 1 it can be easily shown that 
MFPT into the state N x is K ■ (K x + 1) times faster then MFPT into the state N 2 . Here 
K\ = [Ni]/[Ux] is the equilibrium constant between states N\ and U±. When the ground 
state is stable (K ^> 1) the ratio of folding rates into states N x and N 2 is gretaer than in 
the case when ground state is unstable. This is due to the fact that at such conditions the 
state N\ should be considered as a kinetic trap for folding into the state N 2 . 

This conclusion is qualitatively consistent with our numeric results which suggest (Fig. 3) 
that the ratio of MFPT into ground state conformations becomes greater as temperature 
decreases, i.e. both ground states become thermodynamically stable. To test this prediction 
quantitatively, we estimated the temperature dependence of the equilibrium constant Ki 
using histogram technique (Ferrenberg and Swendsen, 1989; Sali et al, 1994; Socci and 
Onuchic, 1995; Abkevich et al, 1995). A long Monte Carlo simulation was performed 
(2 ■ 10 7 steps) in which the conformation shown in Fig. 2a was folded and unfolded many 
times (~ 100) but the conformation shown on Fig. 2b was not yet reached. The statistics of 
occurrence of different states were collected. The two important parameters were taken into 
account: the energy of the chain E and similarity parameter Q\. The quantity of interest 
is the logarithm of the density of states v^E^Qi). Once calculated, it makes it possible to 
calculate population of states at different temperatures. This provided the estimate for the 

3 This is a reasonable assumption because folding rates of sequences, having only one of the 
conformations shown in Figs.l and 2 as native are close to each other. 
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temperature dependence of equilibrium constant K\. Fig. 6 shows that the ratios of MFPT Fig. 6 
into states N x and N 2 at different temperatures fit theoretical curve quite well. 

We have shown that number and strength of local contacts is an important factor which 
determines which of the two global minimum conformations is reached faster. This effect 
was observed for the model with designed and random sequences and for the Go model. 
However the question remains: whether strength and number of local contacts is the only 
factor which determines the ratio of MFPT into ground state conformations. To address 
this issue, we randomly selected two maximally compact conformations with equal number 
of local contacts (six). To ensure that energy of local ground state contacts is the same for 
both structures, we studied their folding using the Go model. We found that folding into 
one these structures is somewhat faster than into another. However, the ratio of MFPT into 
ground state conformations is significantly smaller than what we observed for the structures 
shown on Figs.l and 2 (the data not shown). This ratio never exceeds a factor of three and 
at some temperatures is as small as 30%. Such numbers are typical scatter of MFPT for 
model proteins with different but unique native states. 

This result suggests that for the studied model the number and strength of local contacts 
in the ground state is the only major factor which determines which of the ground state 
conformations is reached faster. 

DISCUSSION 

It is important to compare the properties of natural prions with model polypeptides 
with degenerate ground state which we studied in the present work. Structural studies 
demonstrate that PrP c and PrP Sc differ in conformation (Pan et al, 1993). PrP c was 
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found to have high content of a-helix (42%) and essentially no /9-sheet (3%), whereas PrP Sc 
had a /9-sheet content 43% and an a-helix content 30%. This is consistent with our findings 
because in a /5-sheet non-local contacts dominate, whereas in an a-helix local contacts 
dominate. 

The fact that PrP 5c is much less soluble than PrP c has been well established (Pan et 
al, 1993, Kocisko et al, 1994). As pointed out in the Introduction, one can suggest a "ther- 
modynamic" explanation of prion behavior, namely that the stability of monomeric PrP 5c 
is vanishingly low. The implication of this possibility will be that under no circumstances, 
including any dilution or mildly denaturing conditions can monomeric PrP 5c be observed. 
This seems to contradict the experimental finding that dilution of prion solution from 3M of 
denaturant preserves infectivity while dilution from 6M of denaturant eliminates it (Kocisko 
et al, 1994). However interpretation of these experiments is not entirely clear (e.g. whether 
dilution from 6M of GuHCl leads to folding into PrP c or not) and we cannot completely 
rule out the "thermodynamic" explanation though we find it less likely. Perhaps further ex- 
periments including the stopped flow folding of prions upon rapid dilution from 6M GuHCl 
with subsequent monitoring of secondary structure formation can be useful in clarifying this 
very important question. 

Prions are not the only proteins which can spontaneously undergo global structural 
changes and transform into another stable state. The theory described here predicts that 
in such transitions fraction of local contacts should decrease. What can be observed ex- 
perimentally is the decrease of content of a-helixes accompanied by increase of content of 
/9-sheets. 

It is important to mention that conformational flips from /5-sheet to a-helix are also 
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observed for some proteins (Reed & Kinzel, 1993). However these transitions occur when 
protein environment is altered (change of the solvent or addition of an agent which stabilizes 
a-helical conformation). In this case relative stability of native state can change significantly, 
and the theory presented above is not applicable. We should concentrate on the cases when 
proteins are (metha) stable during their life time in different states under the same conditions. 
For example, upon exposure to bright illumination photosystem 77 irreversibly transforms 
into a stable inactive state. At the same time content of a-helix drops from 67% to 24% 
whereas content of /9-sheet increases from 9% to 41% (He et al, 1991). Another example 
is a human plasminogen activator inhibitor- 1 (PAI-1) which also spontaneously folds into 
a stable inactive state without cleavage (Katagiri et al, 1992). In agreement with the 
presented theory, inhibitory activity of this protein can be restored through denaturation 
and renaturation (Hekman & Loskutoff, 1985; Katagiri et al, 1992). It was suggested that 
transition from active into the latent state is due to transformation of a surface helix into a 
/5-sheet (Mottonen et al, 1992). 

Another interesting question is the evolutionary origin of prions. It was estimated that 
the probability to randomly synthesize a protein sequence with degenerate stable ground 
states by chance is low (Gutin & Shakhnovich, 1993). This may imply that even if some 
properties of the "abnormal" , patogenic, state seem now useless and even harmful, proteins 
could have been specially designed during evolution to have "abnormal" state as well as the 
native one. Finding the possible biological role of the "abnormal" state can be important 
for understanding of protein's properties. Alternatively we can assume that prion proteins 
have only weakly optimized sequences, so that they are intrinsically unstable in their ground 
state conformations at physiological temperature. Then probability of accidental synthesis 

15 



of a protein with degenerate native state is sufficiently high (Gutin & Shakhnovich, 1993). 

Finally we would like to point out to limitations of the present analysis. In order to 
simulate the effect of kinetic partitioning, we had to use a simplified model. Such simplifi- 
cation comes at a price. One limitation is that we studied folding of model proteins much 
shorter than typical experimental systems exhibiting kinetic partitioning. E.g. prions are 
approximately 200 aminoacids long, while in the present study we simulated folding of 27- 
mers and 36-mers. While much longer model proteins can be folded on a lattice (successful 
folding simulations of model proteins of up to 100 aminoacids long have been reported by 
several groups (Kolinski et al, 1993, Shakhnovich, 1994), the detailed analysis presented in 
this paper required thousands of runs to collect sufficient statistics. This is feasible only for 
relatively short model proteins. However, we do not think that this is a crucial limitation 
since comparison of our results with experimental situation suggests that gross structural 
features leading to kinetic partitioning may well be reproduced by simulations of shorter 
chains. 

Another limitation of the present study is that it did not include intermolecular interac- 
tions which are important for transition from PrP c to PrP 5c . Here we should emphasize once 
again that our results are aimed to explain why prion proteins fold into the conformation 
of PrP c after synthesis of their primary structure. Certainly, the presented simulations and 
analysis do not address the mechanism of the conformational change of PrP c into PrP Sc . It 
is possible to simulate small ensemble of lattice chains, or study template-mediated folding 
and address the issue of prion infectivity and conformational transitions (Cohen et al, 1994). 
We are planning to do so in the near future. 
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FIGURE CAPTIONS 



Fig.l (a), (b) Randomly chosen maximally compact 36-mers on a cubic lattice and (c) a 
sequence for which these conformations have equally low energy. 

Fig. 2 Compact 36-mers on a cubic lattice with (a) many local contacts, and (b) without 
local contacts, and (c) a sequence for which these conformations have equally low energy. 
How these structures were found is described in our previous work (Abkevich et ai, 1995). 
Fig. 3 Dependence of MFPT for 36-mers shown on Fig. 2 on the inverse temperature (Go 
model). Gray points correspond to the folding into conformation shown on Fig. 2a, and black 
points correspond to the folding into conformation shown on Fig. 2b. 

Fig. 4 Monte-Carlo folding trajectory at folding transition temperature (T = 0.6) for the 
chain which has two ground state conformations shown in Fig. 3 (Go model), a: the MC-Step 
dependence of the structural similarity Qi to the conformation Aq shown in Fig. 2a. This 
parameter was defined as the number of common contacts in the current and the ground 
state conformation Aq divided by the total number of contacts in Aq. b: The same plot but 
for the conformation N 2 shown in Fig. 2b. 

Fig. 5 Schematic representation of energy landscape in a case of degenerate ground state. 
Aq and A*2 are the ground state conformations, U\ and Ui are unfolded conformations in 
which dominates contacts from the corresponding ground state. 

Fig. 6 Temperature dependence of ratio of MFPT into conformation shown in Fig. 2b 
(MFPT 2 ) to MFPT into conformation shown in Fig.2a (MFPTi). Experimentally observed 
data (grey dots) is taken from Fig. 3. Theoretical curve is calculated with a help of his- 
togram technique (Sali et ai, 1994; Abkevich et ai, 1995). The only fitting parameter was 
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equilibrium constant K which was taken to be equal 81. 
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